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1 Introduction 

Nowadays, the QCD formulated in terms of quark and gluon fields is confirmed as the fun- 
damental theory of strong interactions. However, the quantitative predictions for hadronic 
reactions meet with unsolved problems coming from a non-perturbative character of the 
QCD at large distances. The alternative effective field approach, successfully applied in 
the low-energy hadron physics, is the so-called Chiral Perturbation Theory (ChPT). This 
theory is based on the internal SU(3)l x SU(3)r x U{l)y symmetry of QCD, violated by 
quark masses, and is formulated in terms of hadrons - mesons (as the Goldstone bosons) 
and baryons. The ChPT, more than 30-year history of which can be traced back to the 
pioneering works by Weinberg pQ and Gasser and Leutwyler [2] , has developed into a pow- 
erful tool for investigating the im [3J, ttN [UI5] as well as the few nucleon systems [51 E]- A 
large number of references can be found in the reviews (H |9l [101 [11] . The Lagrangian of the 
effective theory may generally contain infinite number of terms with different number of 
derivatives and unknown constants, and ChPT makes sense only at small particle momenta 
(p) and meson masses (/i), i.e. at p, /i <C A x = 47?/^ ~ 1 GeV, where / ~ 93 MeV is the 
pion decay constant. In this case we have the expansion parameter x ~ /V-A* "C 1, and 
the leading-order (LO) term of the Lagrangian is well denned. 

In this paper we make the first attempt to apply the ChPT to the annihilation processes 
NN — > mesons. We are interested in the case of slow final mesons (with momenta k A x ) 
in the reaction rest frame. This situation takes place, for example, in the multi-pion 
annihilation process NN — )• mr with large n ~ y/s/fx ^> 1. The appropriate data up to 
n = 9, obtained from pp annihilation at rest can be found in Ref. [12]. Theoretical study 
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of such processes involves calculation of a large number of diagrams and integration over 
the final many-particle phase space. In this respect, the reactions with a small number of 
heavier final mesons r\ and K(K) seem to be more attractive, and corresponding data are 
available as well [13J. In this paper, we consider the pp annihilation processes 

(a) pp — > TjTjTj, (b) pp — > r]K°K°, (c) pp — > r]K + K~ (1) 

at rest. Here, the typical meson momentum is k ~ 300MeV/c, the ChPT expansion pa- 
rameter value x ~ /V^x ~ 0-5, and we may expect a sensible results from LO calculations 
in the effective chiral theory. In the present paper, we perform calculations for the proba- 
bilities (branchings) of the reactions ([I]) at rest in the tree approximation, using the SU(3) 
chiral effective Lagrangian at LO [91 [TO]. 

The paper is organized as follows. In Section 2, we describe the SU(3) chiral Lagrangian 
and write down its parts used in the tree calculations. In Section 3, we calculate the 
amplitudes for the reactions ([JJ). In Section 4, we present our numerical results for the 
branchings, compare them with the experiment and discuss the results. Section 5 is the 
Conclusion. 



2 Lagrangian 

We consider the SU(3) chiral effective Lagrangian at leading order. The meson part £ 
of the Lagrangian can be written in the form 




Here: tt is the octet matrix of the meson fields, and U represents their customary ex- 
ponential parametrization; A a (a = 1, ...8) are the Gell-Mann matrices, normalized as 
(AaA;,) = 25 a b', the angle brackets (• • •) denote the trace over flavours; / ~ 93 MeV is the 
pion decay constant; M = diag(m u , m d , m s ) is the quark- mass matrix. The first term C\ in 
Eq. d2J) is the chiral symmetrical part of the Lagrangian C nn . The second term £ 2 contains 
the quark-mass matrix M and spontaneously violates the chiral symmetry. The parameter 
u in the term £ 2 relates the meson and quark masses in the form 

ml = a(m u +m d ), m 2 K ± = a{m u + m s ), 

m KO,Ro = &( m d+ m s ), m\ = \a(m u + m d + 4m s ), 
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known as the Gell-Mann, Oakes and Renner relations [H].0 To obtain Eqs. ([3]) one can 
expand the Lagrangian £ 2 © in the meson fields and identify the quadratic term for a 
given meson 7c a with the mass term — }^m 2 'n 2 a . 

The baryon-meson part of the effective Lagrangian can be written as 



£-kB — £3 + ^4) 

C 3 = (B(iVrf-m B )B), C i = (BY l5 {D{A^B} + F[A^B])). 



(4) 



where 



V ll = d ll B + [V li ,B], 



A, 



B = -=J2 B a( x ) X ° 



V2 



s- 



27T 

S" 1 



(5) 



17, 



P 



n 

v ~ x/rv 

Here: S is the octet matrix of the baryon fields; m B is the baryon mass; D and F are the 
usual axial-vector meson-baryon coupling constants; [A, B] and {A, B} mean commutator 
and anticommutator for the operators A and B. 

We expand the Lagrangians (j2J) and (jlj) in meson fields and retain only the 
interaction terms of interest for the pp — >(3 mesons) amplitude in the tree approximation. 
So, we use the terms of the fourth order in meson fields from £ T7r and the terms up to the 
third order from C^b- Then, the total interaction Lagrangian takes the form 



Li 



L = L x + L 2 + L 3 + L 4 , 



(6) 



3J2 N L ? "J LA* ' J/ ? 



L 3 = i(%^, J B]), 
1 



2f 



L 4 = £4 (see Eq. g 
A, 



4^7T- — [[^,77],^], 



Finally, since m u ^ <C m s , we take the u and d quarks to be massless. Thus, we have 
M = diag(0, 0, m s ) and o = 3m^/4m s (see Eqs. (J3])) for the term L 2 in Eqs. (jUJ). 



1 These also include the equation for the quark condensate (qq) = — 3a f 2 
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3 The amplitudes 



The full sets of the tree diagrams for the reactions (Ilk), and ([Tb,c) are shown in Figs. [T] 
and |2l respectively. Here we write down the corresponding amplitudes. The vertices of 
these diagrams with a given mesons and baryons are fixed by the Lagrangians L]_4, and 
are given in Appendix A. We shall use the notations: 

jf = tf% Gb(p)= * + ?* , u=V*ii(*), v=V2m~(0,- X + ). (7) 

Hereafter: Gb{p) is the propagator of baryon with 4-momentum p; u(v) is the bispinor of 
the initial proton(antiproton) at rest; (p and \ are the spinors ((p + (p = x + X = 1); m is t ne 
proton mass. 

3.1. The reaction pp — > rjrjr] 

The amplitudes for the reaction (JT^i) read 

%M A = V^vft 3 7 5 G p (k 3 -p 2 )lik2j5Gp{Pi-kiWi r y5U + (k^h permutations), 



PPV 

iM B = V m V^XG v , X= vP l5 u = -Am 2 (xM- 



Here: hi = (oji, fcj) are the 4-momenta of the final r] mesons; p±(p2) = {m, 0) is the 4- 
momentum of the initial proton(antiproton) at rest, and P = p\ + P2 is the total 4- 
momentum; G p (- ■ •) is the proton propagator (5 = p, see Eq.(J7J)); G v = l/(4m 2 — m^) 

is the 77 propagator in the diagram (B); and V\ are the vertex constants, given in 
Appendix A. The amplitude Ma flU) contains 3! = 6 terms, arising from permutations of 
77's momenta fcj. Further, we express the amplitudes M (|SJ) in the form 

zM = + io-B)^, (9) 

making use of Eqs. (JTj). Then for the diagrams in Fig. [I] we obtain 

(A) : A,B = (see Appendix B); 

(2) ( 10 ) 

(B) : A = -Am 2 V m V^G v , B = 



We also consider the simplified amplitudes, calculated with final mesons at rest. In this 
"threshold" s-wave approximation we take 

T 

ki = 0, Ui = fii + —, T = 2m - iii - /i 2 - /i 3 , (11) 

where //j are the meson masses, and T is the total kinetic energy (the final mesons are 
produced at rest as if their masses are /2j = ui). In this case we have two graphs 

V V 3 - 2mii) 3 

(A) : A<°> = — ^ rf; (B) : A<°> = A(Eq. (HDD). (12) 

(2m — a;,,)" 1 
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and the terms of the type icrB vanish (we use a superscript "(0)" for the amplitudes in 
the "threshold" approximation). 

3.2. The reactions pp — > r)KK 

The diagrams for these reactions are shown in Fig. [2j Let us denote here the 4-momenta 
of final rj, K and K by k\, k 2 and k^, respectively. Making use of the Lagrangians, given 
in Appendix A, we obtain the pp — > rjKK amplitudes in the form 

iM Al = V ppv V 2 Ak v^ b GA{h-p2)hl^ G p{Pi- k i)hl^ u ^ 

iM M = V ppv V 2 Ak v^ 1 ^G p (ki-p 2 )%^GA{pi-k 2 )^2lbU, (13) 

%M M = V AAv V 2 AK v^ 5 G A (k 3 -p 2 )^ 5 G A (pi-k 2 )^5U] 

iM B (l, a) = 2V la (2mu v + m 2 K - m 2 KR )XG a , 

iM B (2, a) = V 2a XG a , (a = V , tt°), V ia = V ppa V^= R 2) ; (M) 
iM c = -V c v(2# 1 -fa-fo)'y5U, V c = V mKR ] (15) 



(16) 



iM Dl = -V D vk , 3 j 5 G A (k 3 -p 2 )(tk 2 -k'i)u, 

iM D2 = -V D v{^ 1 -%)G A {p 1 -k 2 )^2lbU , V D = V pAK V pAKr) 

iM Fl = -V F v^ b G A (ki-p 2 )(^ 2 -^)u, 

iM F2 = -V F v{ty 2 -k' z )G A {p l -k 1 )k l 1 ^u , V F = V m V ppKR . 

Here: m KR is the effective mass of KK system; G a = (4m 2 — m 2 )' 1 is the meson propagator 
(a = rj, 7T°), and m a is the meson mass. The amplitudes M B (i, a) ( I14p are specified by the 
Lagrangian Li (i — 1,2), determining the 4- meson vertex of the diagram (B) in Fig. |2j and 
by the exchanged meson a — rj, 7r°. Representing the amplitudes in the form fl9]), we obtain 



(^1,2,3, D 1>2 , F 1;2 ) : A,B = (see Appendix B); 
(S li0 ) : A = -8V la m 2 (2mu v + m 2 K - m 2 KR )G a ; 
{B 2A ): A = -4V 2a m 2 G a , 

(C): A = AV c m 2 (2oo v -uj K -uj R y, B = for (B), (C); 
In the " threshold" approximation ( II ip we have B = and 

(A )■ A (0) - V ppTi V pAK 2muj yU 2 K (Q) _ V AAv V 2 AK 2muj ri uj 2 K _ 

1 1,2) ' (2m-m ri )(M A +m-m K y 1 3j ' (M A +m-m K ) 2 ' 

A^=-8V la m 2 (2mio v -3uj 2 K )G a , A<°>= (see Eq. flE}); 

(C): AM=8V c m 2 (u v -u K ); (D^): A (o) = y ^i^-u,) _ {F):A w =Q 

M A +m — 0OK 



(17) 
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3.3. A remark on power counting 



Let us look at the tree approximation in connection with the ChPT power-counting rules 
and give the dimensional estimation of the amplitudes. We express the amplitudes through 
the following factors: V (the product of the vertex constants); 2m (normalization factor 
for bispinors, i.e. uu = 2m); u n , where u is typical value of the meson total energy (for 
the meson-baryon vertex of the diagrams (B) we replace u — > 2m) and n is the number of 
derivatives in the Lagrangian, determining a given vertex); 1/u and 1/Am 2 for baryon and 
meson propagators, respectively. Let us rewrite the factors V through the dimensionless 
constants V. Then, V = Vm 2 /f 3 ~ Vu 2 / f 3 for diagrams (7> 2i(l ) in Fig. [2] (and diagram 
(B) in Fig. [TJ, and V = V / f 3 for all the rest of diagrams. The resulting estimations of the 
diagrams in Figs. [T] and |5] are 

(A) ~ (C) ~ (D) ~ (F) ~ 2mu, (B) ~ u 2 . (19) 

Since u ~ 2m/ 3 ~ m, the amplitudes ffl9|) are approximately of the same order of mag- 
nitude. However, the V values, fixed by the SU(3) Lagrangian, are very different (these 
numbers are presented in Table 2). Thus, the contributions of the tree diagrams are very 
different, especially in the pp — > r]K + K~ CcLSC, clS discussed below. 



4 The results 

Applying the amplitudes of sect. 3, defined by Eq. (EJ), we get the pp-annihilation cross 
section va at rest {y — > 0) in the form 

va( P p^...) = —J\M\ 2 dT 3 , d7S = -^_ ( y T s = 2m). (20) 

Here: / is the identity factor, and 7 = 1/3! (7 = 1) for the 777777 (j]KK) channel; \M\ 2 = 
\{\A\ 2 + \B\ 2 ) is the amplitude flU]), squared and averaged over the initial p and p spin 
states; dr 3 is the phase-space element (in the Feynman normalization) of the final 3-meson 
state, written for the unpolarized case; k\ = Q(y/s, fix, w) and q = Q(w, /j 2 , ^3), where 
Q(m,x,y) is the relative 3-momentum in the pair of particles with masses x and y, and m 
is the effective mass of this pair; w is the effective mass of the final meson pair (2+3), and 
z = cos 8, where 9 is the polar angle of the relative motion in the (2+3) pair. The results 
obtained from Eqs. ( 120]) with the amplitudes, given by Eqs. ( |T0|) or ( TT71) . will be referred 
to as the "exact" results. 

We also give simplified ("threshold") version of calculations, making use of the ampli- 
tudes, given by Eqs. (fl2~|) and ([18]) . Then, the cross section can be written as 

877-r 6471-^ (/ii +^2+^3) 7 
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where is the sum of terms, given by Eqs. (fT2l or (TIB"]) for 777777 or r]KK channels, 
respectively; is the non-relativistic phase space [15] of final 3-meson state, and T = 
2m — fix— /x 2 — /X3 is the released energy. 

Further we give the results in terms of branchings (Br), where 

Br(pp-+---) = va &P^'") vaZ 1 (v ^ 0) = 36.5 mb. (22) 

PP 

The total pp-annihilation cross section va^ n can be obtained from the imaginary part of 
s-wave scattering length Ima s = —0.69 fm, extracted in Ref. [IE]. Then, we get va^ n {v — > 
0) = (8n/m) Im(-a s ) = the value in Eqs. (|22j) . 

Our numerical values of branchings as well as the experimental ones are presented in 
Table 1. We give both "exact" and "threshold" results (the latter are in brackets) for two 
sets (a)- (b) of values of D and F constants, which enter the term L4 of the Lagrangian (J6l) 
and were determined from semileptonic baryon decays. These sets are 

(a) D = 0.76, F = 0.48 [13; 

(23) 

(6) £> = 0.80, F = 0.46 [IB]. 

The results in Table 1 show some sensitivity to D and F values. Let us consider different 
channels separately. 



Table 1: The branchings for pp-annihilation reactions at rest. The data [13] 
were obtained in liquid H 2 target. The data set for 3?7 channel was extracted 
from pp — > 67 and pp — > IO7 events. The K°K°r] branching is obtained from 
K s K s r] data [33] as Br(K°K°r]) = 2 Br(K s K S T]). 



Channel 


Br (theor)-lO 4 


Br (exp)-10 4 [13] 


pp — > r\r\r\ 


(a) 0.424 (0.343) 
(6) 0.314 (0.275) 


4.22 ± 0.24 ( 67) 
3.93 ± 0.68 (IO7) 


pp — > r]K°K° 


(a) 1.488 (0.763) 
(6) 1.676 (0.510) 


2- (2.2 ±0.7) 


pp — > 7]K + K~ 


(a) 182. (0.0035) 
(6) 157. (0.0011) 





pp — y rjrjr). The calculated branchings underestimate the data [13] approximately 
by one order of magnitude. Let us consider our result (a) in more detail. The "exact" 
version of branching value exceeds the "threshold" one by ~ 24%. The diagrams (A) and 
(B) of Fig. [U taken separately, give the "exact" ("threshold") Brirjrjrj) values 4.65 • 10~ 8 
(1.16 • 10~ 6 ) and 4.50 • 10~ 5 (4.81 • 10~ 5 ), respectively. Thus, the dia gram (B) dominates 
and (A)-(B) interference is destructive. The double-baryon-exchange diagram (A) give 



7 



very different and relatively small contributions in the "exact" /"threshold" versions. The 
amplitude Mg is approximately the same in both versions^] 

pp — >■ rjK°K°. Our results (a)-(b) underestimate data by a factor ~ 3. The ex- 
perimental branching should also include the channels with direct resonance production. 
One of them is the 4>r) channel, studied in the reaction pp — > i]K + K~, and the extracted 
branching Br( 3 S u pp_^ <fa K + K~ri) = (0.76 ± 0.31) • 10" 4 pH EU [21] , which should be 
the same in the K K°rj channel. Substracting this value with weight factor |, we obtain 
(2 ■ 2.2 — | 0.76 = 3.83) ■ 10~ 4 , what is closer to our results. Some contribution should 
also come from the / Q (1500)?7 channel [22] due to the decay / Q (1500) ->• KK [23]. With 
these reservations we may accept the approach based on the diagrams of Fig. [2] as the 
"background" model for the pp — > t]K°K° reaction at rest. 

Contributions of the diagrams in Fig. [2] to the branching as well as signs of the cross 
terms are very different. The largest contributions come from the diagrams (_B li7r o), {D\ + 
D 2 ) and (Fi+F 2 ), and corresponding Br(K°K°r]) ■ 10 4 values for variant (a) are 2.30 (0.39), 
3.09 (0.016) and 2.45 (0), respectively. The results of "exact" and "threshold" versions are 
drastically different, and the (-^1,2) terms vanish in the latter case. 

pp — > rjK + K~ . For this reaction the calculated branching exceeds the result for 
Br(K K°i]) approximately by 2 orders of magnitude. The experimental value for the 
charged channel case is unknown. We roughly estimated it from the data on the annihila- 
tion frequency Y [20]. For the s-wave branching we may write Y > (1 — f p ) ■ Br(K + K~r]), 
where Br(K + K~7j) = |Br( 1 S'o) + |-Br( 3 Si), and f p is the fraction of annihilation from 
the p-wave statesfl Taking the values Y = (8.17 ± 0.37) • 10~ 4 [20] and f p = 0.13 [H] for 
pp — > f]K + K~ at rest in the liquid target, we obtain an estimation Br(K + K~i]) < 9.4- 10~ 4 
for the upper value. Thus, the theoretical value exceeds the estimated data by more than 
one order of magnitude. 

It is worth noting that the "exact" version of our results exceeds the "threshold" one 
approximately by 5 orders of magnitude, i.e. the latter approximation does not work at 
all. Very small Br(K + K~rf) values in the "threshold" version are due to cancellation of 
different terms in the total amplitude. Table 2 contains the branching values, obtained 
from different diagrams of Fig. [2] separately. Here, we give the "exact" and "threshold" 
results (the latter in brackets) for variant (a) and absolute values \ V\ of coupling factors V 
(see Eq. f[T9~|) ) for each diagram. The subscripts S° and A specifies the exchanged baryon A 
in Fig. [2j Further, we discuss the results of the "exact" version. 

The largest contributions come from the diagrams (Ax+2,a), (^3, a), (-Di+2,a) and 
(F 1+2 ). Note that the terms (-Di+2,a) have the largest factor \V\ = .238 due to the large 
coupling constants, containing A baryon. With these terms excluded, i.e. for the incomplete 
set (^4 1+2+ 3 s°) + C^) + (C) + (D 1+2 s o) of diagrams, we have Br(K + K~rf) • 10 4 = 6.2. Suc- 
cessive addition of the terms (Fi +2 ), (^1+2+3, a) and (-Di+2,a) give the values Br(K + K~rf) ■ 

2 The difference of Br(rjrjrj) values from diagram (B) comes from the phase space difference 73/73°' = 
4.50/4.81, where r 3 (0) is the non-relativistic phase space in Eqs. (pH]) . 

3 Here, we neglect the correction factors E^Sq) and E( 3 Si) for the statistical weights \ and | of 1 Sq 
and 3 Si states, equating them to 1. For example, E( 3 Si) = 0.989 for the liquid H 2 target [T9"] . 



8 



10 4 = 12.1, 42 and 182, respectively. Note that the cross terms from these contributions 
are constructive. Thus, we obtain a large branching value Br(K + K^rf) ^> Br(K°K°r]), 
and the main reason comes from the diagrams containing A exchange, which are absent in 
the neutral K°K°r] channel. 



Table 2: The branchings for pp — > r)K + K at rest from the diagrams 
of Fig. [2j The subscripts S° and A specifies the exchanged baryon A. 



diagram 


Br ■ 10 4 


\v\ 


diagram 


Br ■ 10 4 


1^1 


0*4-1+2, £°) 


.0092 (.0029) 


.0038 


(4»,E°) 


.0095 (.0023) 


.0086 


(A+2, A) 


4.57 (1.35) 


.079 


(A3, a) 


5.52 (1.19) 


.177 




0.91 (0.16) 


.049 


(Bi, w o) 


2.55 (0.46) 


.089 




0.44 (0.46) 


.037 


(B 2>7r o) 


0.14 (0.14) 


.045 


(C) 


0.44 (0.15) 


.069 


(A.+2,s°) 


0.83 (.0047) 


.030 


(A + 2,a) 


61.9 (0.32) 


.238 


(F 1+2 ) 


10.6 (0) 


.098 



5 Conclusion 

The pp-annihilation processes pp — > r\r\r\ and pp — > r]KK at rest were considered in the 
tree approximation in the framework of the SU(3) chiral effective Lagrangian at leading 
order. The following results were obtained. 

1) The calculated branchings for the r\r\r\ and K°K°7] channels are several times smaller 
than the experimental values. Our results underestimate the data approximately by one 
order of magnitude in the 777777 case and - by ~ 3 times in the K K i] one. One of the 
important missing contributions may be the direct production of intermediate heavy reso- 
nances, not incorporated into the pseudoscalar SU(3) octet of mesons. Thus, for example, 
an essential role of the resonance /o(1500) in the 7777-mass spectrum was found in the reac- 
tion pp — y rjrjr) at beam momenta 600 — 2410 MeV/c [22J. Due to a large /o(1500) width 
T w 109 MeV [23], one may also expect this resonance to affect visibly the calculated cross 
section of the reaction at rest. The other intermediate heavy resonances can also con- 
tribute. On the other hand, we can not exclude a sizeable role of the next-to-leading-order 
(NLO) corrections, including the loop diagrams. Potential pp-interaction in the initial state 
can also change the reaction cross sections calculated in the ChPT approximation [see, for 
example, Ref. [24]]. 

Thus, we consider our estimation of the tree diagrams as the calculation in the Born 
approximation and as the background model of the given processes. 

2) The calculated branching for the K + K~rj channel exceeds the estimated data by more 
than one order of magnitude. Formally, the data on this branching are unknown, and our 
statement is based on rough estimation of the "experimental" value, given in Sect. 4. 
Perhaps, the possible resonance-production mechanisms as well as the NLO corrections 
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could partly explain this strong discrepancy. However, here the results obtained seem to 
be against the applicability of the chiral effective approach to the annihilation processes. 
An important reason of large calculated K + K~rj branching comes from the large SU(3) 
vertex constants, involving A-baryon exchange (not occuring in the neutral K K i] case). 
Perhaps, the NLO corrections may be also enhanced due the A vertices and improove the 
theoritical branching value. This is an opened question, faced to the further investigations. 
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Appendix A 

Here we write down the terms of the Lagrangians Li_ 4 fl6]), corresponding to the vertices 
of the diagrams in Figs. ([T]) and (j2J) with a given baryons (B) and mesons (tt). 

1. The four-meson terms of Li (JH]): 

Li(r]T] KK) = V [l) KR [r]d^ (Kd»K+ Kd"K) - d^d^KK- ^ d^d^K], 



r]r]KK 

Li(-K°r]KK) = V^ vKR [(-K% V + V dy)(Kd^K+ Kd»K) 

- 2<V/ d^°KK- 27 1 TT°d IM Kd^K], 

V {1) -— v {1) =-v {1) - 1 

v riqKK 4 p ' TT°r)K+K~ V w a r)K°K ^y/Sf 2 ' 



(A.l) 



2. The four-meson terms of L 2 ©: 
L 2 (n a n b 7i c n d ) = — ^f|^ d 7r a 7r 6 7r c 7r d , 



nl 



T/( 2 ) - 41 ^L. T/(2) - 9! _ _t/(2) _ m \ 



(A.2) 



Here, we introduce the identity factor nl, where n = 4(^ 4 ), 2(j]T}KK), l(7T°r]KK). 
3. The terms of L3 ([6]) for B a B b — > 7T c 7id'- 

L 3 (B a B b ir c iT d ) = iV Ba§b7Tc7Td u&7/i"a Od<9 M 7r c - 7r c <9 M 7r d ), 



V v ^ v 1 v 1 

V pTPK+rj — ~ V Z°pK-r) — -^pi Vpp K + K —J) v ppK°K° — J7j 

4. The terms of L 4 © for B a B b ir c : 

Li(B a B b ix c ) = V BaBb7Tc u blt a h u a d"Tt a , V BaBb7Tc = V BbBa - c ; 
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(A.3) 



v - 3F ~ D v - D + F v- - D ~ F 

v PPn ~ 2v /3j ' pp*" ~~ 2f ' v p x°k+ - 2 j ' 

D — F D + 3F D 

Vpt]+K0 ~ ~w pRK+ ~ ~^W' v ^ v = ~ VaKv = Wf 



(A.4) 



5. The terms of L 4 ([6]) for pp — )■ rjKK: 

U{pp V KK) = V mKR u pl ^u p [2KKd^ - V (Kd»K + Kd^K)}, (A.5) 

Up{up) is the Dirac spinor of the proton (antiproton) in items 3—5. 

Appendix B 

Here we give the expressions for the amplitudes M^(jSj), M^ 123 (ITBl . Md 13 (ITS]) and 
M Fl 2 (ITS]) , represented in the form (jUJ) . 

1. The amplitudes and M^ 123 . For the amplitude, say M Al , we obtain 

(Ax): A = c[a + b(k 1 k 3 )}, B = cb[k 3 xk 1 ], (B.l) 



where 



c = 2mV D^D 3 \ D li3 = m 2 - 2mwi, 3 + /i 2 3 - M 2 3 , 



a = B 3 (A 1 kl + 5iw 2 ) + 5iA 3 fc 2 , b = A 1 (B 3 + A 3 tu 2 ) + B 1 A 3 , 

A,3 = m + M 1)3 , Si i3 = (Mi j3 - m)wi i 3+ M?, 3 - 

Here: V = VppnVpAK * s ^ ne ver tex production; k i: uji and /ij are, respectively, the 3- 
momentum (&j = |fcj|), total energy and mass of the final meson (see Fig. [2]); M 13 are the 
exchanged-baryon masses (Mi = m and M 3 = for the diagram A\ on Fig. [2]). The 
amplitudes M A (Fig. H]) and the rest of the amplitudes M A . (Fig. [2]) can be obtained from 
Eq. (B.l) by the proper replacement of the vertex factor V, final-meson momenta and 
baryon masses Mj. 



2. The amplitudes Md 12 and Mp l a . For these amplitudes we obtain 
pi, 2 ): A = c li2 [a li2 +6(A; 2 -A; 3 2 )], B = 2c 1)2 6 [fc 3 x fc 2 ] ; 

(F 1)2 ): J 4 = c[(w 3 ,2-W2>? + 2m(fc 3 V fc 2 2 i]. B = 4cm [fc lj3 x fc 2 ,i] 

Here: 

c 1)2 = -2mF D Z}^, D lj2 = m 2 - 2mw 3j2 + /i 3 2 - m A , 
ai,2 = (wi-w 2i3 )[(m j4 -m)a;3 i2 +/i 2 i2 ], & = m+ m A , 
c = — 2mVp(/i 2 — 2mwi) _1 ; 
Vd and Vf are the vertex productions (IT6]) . 



(B.2) 
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Figure 1: Feynman tree diagrams for the pp — » amplitude. Solid and dashed lines 
correspond to the baryons and mesons, respectively. 
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Figure 2: Feynman tree diagrams for the pp — » rjKK amplitude. Solid and dashed lines 
are the same as in Fig [TJ The exchanged baryon A is S + (S or A) in the case of the 
7]K°K° (r]K+K-) channel. 
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